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20.  ABSTRACT  (Canttmm  on  ravora*  N*  If  mtMMl)  aM  lAnKfy  by  Mac*  anNr) 

Coherent  microwave  backscatler  is  the  basis  of  a probe  technique  useful  for  studying  surf  ce 
water  waves  in  a laboratory  tank.  A parabolic  antenna  is  focussed  to  give  a plane  wave  at  short 
range,  typically  1.2  m,  and  an  illuminated  area  of  control'able  size.  The  technique  strongly  dis- 
criminates against  all  water  waves  except  those  which  have  a particular  Bragg  resonant  wave'ength 
ranging  between  0.25  cm  and  10  cm  and  propagate  parallel  to  the  plane  of  incidence.  The  resulting 
dopplcr  spectra  provide  a powerful  tool  for  wind  wave,  breaking,  and  spray  studies.  A unique 
method  'or  directly  measuring  the  water  wavenumber  resolution  is  described,  together  with  a 
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technique  for  absolute  calibration  of  the  backscattered  power  in  terms  of  wave  height.  As  a con- 
sistency check  cm  all  the  calibration  procedures,  the  angular  dependence  and  absolute  value  of  the 
first-order  scattering  cross  section  from  monochromatic  water  waves  agrees  with  prediction. 
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WIND  WAVE  STUDIES:  PART  2— THE  PARABOLIC  ANTENNA 
AS  A WAVE  PROBE 


INTRODUCTION 

A wa  ;e  probe  which  uses  microwave  doppier  backscatter  has  been  developed  and 
used  for  several  years.  The  description  given  here  k intended  to  make  this  technique 
for  studying  water  waves  more  widely  known  and  more  accessible  to  other  groups.  The 
general  technique  of  doppier  backscatter  is  discussea  in  section  1,  as  ate  several  unique 
advantages  of  the  wave  probe.  Section  2 treats  the  illumination  pattern  from  a focussed 
parabolic  antenna,  and  section  3 treats  the  backscatter  response  of  such  an  antenna.  A 
first-order  derivation  of  the  scattering  produced  by  a rough  surface  is  described  in  the 
appendix.  These  treatments  are  intended  to  show  how  the  illuminated  area  and  the 
water  wavenumber  resolution  can  be  selected  by  appropriate  choice  of  antenna  proper- 
ties. A unique  method  of  directly  measuring  wavenumber  resolution  is  described  in 
section  4;  representative  results  are  presented.  Measured  real-space  antenna  patterns 
give  the  illuminated  area  at  the  water  surface,  which  is  necessary  for  absolute  calibration 
of  the  backscattered  power  in  terms  of  wave  height.  Backscatter  produced  by  a known 
monochromatic  wave  is  compared  absolutely  with  the  return  predicted  by  the  first-order 
theory  as  a consistency  check  on  all  the  calibration  procedures  described. 


1.  WAVE  PROBE 

Waves  on  an  air/water  interface  can  be  directly  studied  using  the  microwave  doppier 
backscatter  technique  because  of  the  large  microwave  reflection  coefficient  at  this  inter- 
face [1] . Microwave  radiation  is  focussed  on  the  surface  using  a parabolic  antenna  at 
short  range,  where  “short”  will  be  defined  later  in  this  section.  The  radiation  scattered 
directly  back  into  the  illuminating  antenna,  or  backscatter,  is  carefully  measured  for  both 
amplitude  and  doppier  frequency  shift.  In  other  words,  the  backscattered  signal  is  coher- 
ently detected  using  a portion  of  the  transmitted  signal  as  a local  oscillator.  Hie  scatter- 
ing geometry  is  illustrated  in  Fig.  1.  Microwave  scattering  from  small  amplitude  Water 
waves  is  discussed  in  detail  in  the  appendix. 

For  incidence  angles  0 away  from  0'  and  90°,  essentially  all  of  the  first-order  back- 
scattered  power  is  produced  by  water  waves  which  are  traveling  in  the  direction  defined 
by  the  (Jane  of  incidence  and  which  have  a resonant  wavelength  given  by  the  microwave 
wavelength  Xo  and  8 [2].  More  specifically,  the  waves  contributing  to  first-order  back- 
scatter must  satisfy  the  two-dimensional  Bragg  condition 


Note:  Manuscript  submitted  November  27,  1974. 
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kx  - 2kQ  cos  fl  = kB 


(1) 


and 


ky  = 0 (2) 

A wavenumber  k has  the  usual  definition  2*7>.q,  that  is,  kQ  is  given  by  2ir/X^.  The  x 
direction  lies  in  the  plane  of  incidence,  along  the  undisturbed  surface,  as  indicated  in 
Pig.  1.  Equation  (1)  reduces  to 


, - ^0 
B " 2 cos  d 


(3) 


Ihe  region  in  k space  singled  out  for  study  is  illustrated  in  Fig.  2,  where  resolution  is 
indicated  by  and 


Ibis  Jfe-space  selectivity  is  possible  because  the  measurement  is  inherently  area 
extensive.  In  a sense,  one  is  using  an  infinite  number  of  point  probes  spread  in  a dense 
array  over  the  illuminated  area  Hie  field  distribution  over  this  area  at  the  water  surface, 
and  hence  the  resolution  about  k = (kg,  0),  is  determined  by  antenna  geometry.  Differ- 
ent values  of  ks  can  be  chosen  by  varying  8 or  Xo  or  both.  Our  measurements  extend 
from  a water  wavelength  of  10  cm  (4.30  GHz,  70°)  to  0.25  cm  (70.1  GHz,  30°).  Spatial 
averaging  over  a section  of  tank  width,  common  to  point-probe  techniques,  is  done 
automatically. 
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Fig  2— Region  in  water  wave  k apace  about  (kg,  0) 
studied  by  this  probe  technique.  Resolution  A kx,  Aky 
is  discussed  in  the  text. 


The  output  of  this  wave  probe,  when  run  through  an  audio  spectrum  analyzer,  gives 
the  wave  spectrum  evaluated  at  (kg,  0),  0,  to).  Here  to  is  2*  times  the  water  wave 

frequency,  and  , ky , w)  is  the  Fourier  transform  of  the  surface-displacement  covari- 
ance, as  defined  by  Phillips  [3].  This  output  is  an  instantaneous  doppler  spectrum  and 
is  usually  averaged  in  time  to  give  mean  doppler  spectra.  The  region  in  water  wave 
(cd,  kx)  space  probed  is  shown  in  Fig.  2b.  This  spectral  information  is  available  because 
the  measurement  is  temporally  extensive.  Spectral  resolution  is  limited  by  data  acquisi- 
tion time  and  ultimately  by  the  antenna  pattern. 


To  further  clarify  the  concept  of  a doppler  spectrum,  consider  the  doppler  shift  pro- 
duced by  a mechanically  generated,  small  amplitude,  monochromatic  wave.  Refer  to 
Fig.  1,  with  such  a wave  propagating  through  the  illuminated  area.  Each  point  on  the 
waveform  surface  has  a velocity  component  vd  along  the  microwave  line  of  sight: 

vd  = cos  0 , (4) 

where  up  is  the  water  wave  phase  velocity  in  laboratory  coordinates.  Note  that  the  scat- 
tering occurs  at  the  air/water  interfx-e,  or  waveform  surface,  and  not  from  individual 
water  “particles”  or  from  some  “packet”  of  waves.  Thus  in  Eq,  (4)  the  phase  velocity 
is  the  correct  quantity,  not  a particle  orbital  velocity  or  a wave  group  velocity.  The 
microwave  doppler  frequency  shift  A f relative  to  the  incident  microwave  frequency  is 


Af  __  2ud 
fo  c 


where  c is  the  speed  of  light.  Using 


w * vpk, 

which  is  true  for  all  wate:’  waves,  in  Eq.  (5)  gives 


(5) 


(6) 
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® f water'  <7> 

The  doppler  frequency  shift  of  a first-order  Bragg  resonant  water  wave  is  always  the  fre- 
quency of  tike  water  wave  in  the  laboratory  frame. 

As  part  of  the  antenna  response  measurements  to  be  discussed  Later,  the  doppler 
shift  produced  by  a mechanically  generated  monochromatic  wave  was  measured  directly 
for  several  values  of  kg  The  results  are  given  in  Fig.  3.  For  these  data,  Xo  was  held 
constant  at  3.20  cm  (9.375  GHz),  and  kg  was  changed  by  varying  6 in  5°  increments. 

In  each  case  the  water  wave  frequency  giving  maximum  backscatter  is  plotted  vs  the  fre- 
quency determined  from  kg  using  the  dispersion  relation  for  infinitesimal  water  waves: 

w2  = ^980  j k + ^75-^-^  fc3  (81 

As  seen  in  Fig.  3,  agreement  is  good.  The  slight  systematic  deviation  may  be  due  to  an 
inaccurate  0 calibration  or  to  borcsight  error. 


CALCULATED  BRAGG  RESONANT  FREQUENCY  (Hr! 


Fig.  3 — Experimental  verification  of  the  fact  that  the 
dopple.-  shift  of  a monochromatic  Bragg  resonant  wave  is 
given  by  the  wave  frequency.  The  measured  doppler 
frequency  ia  plotted  vs  the  water  wave  frequency  calcu- 
lated from  the  known  value  of  * g 


Several  other  features  of  this  wave  probe  should  be  mentioned.  Dynamic  range  is 
limited  by  the  audio-frequency  signal-processing  equipment  A range  of  103  in  wave 
hei^tt,  or  60  dB  in  backscattered  power,  is  readily  achievable.  Sensitivity  is  limited,  in 
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practice,  by  how  well  the  unshifted  return  from  stationary  room  object*  can  be  nulled  on 
a microwave  bridge.  Surface  displacements  of  10~2  mm  are  readily  detectable,  and  line- 
arity becomes  a problem  only  for  large-amplitude  waves.  Microwaves  are  “contactless.” 

A slight  beating  is  the  only  perturbation  of  the  water  surface  produced,  and  this  is  usually 
negligible.  Thus  small  waves  can  be  accurately  studied  in  the  presence  of  longer  and  larger 
amplitude  waves.  Because  water  droplets  scatter  microwaves  via  a different  mechanism, 
wave  breaking  and  the  presence  of  spray  can  easily  L-e  detected  (4) . Representative  dop- 
pler  spectra  from  wind-generated  waves  are  shown  in  Fig.  4.  The  microwave  wavelength  is 
7.00  cm  and  d is  30°.  The  air  friction  velocity  p*  varies  from  7 to  124  cm h correspond- 
ing to  wind  speeds  from  1.5  to  15  m/s.  Interpretation  of  the  complex  features  appearing 
at  higher  winds  is  beyond  the  scope  of  this  report.  Known  sources  of  doppier  bandwidth 
are  finite  scatterer  lifetime,  advection  of  scatterers  by  larger  waves,  and  higher  order 
scattering. 
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Fig  4 — Representative  doppier  spectra  fcr  several  wirrd  speeds.  lit  each  esse, 

f0  = 4.3  GHi  and  9 = 30° 


2.  DERIVATION  OF  THE  FIELDS  AT  THE  WATER  SURFACE 

Because  of  restrictions  due  to  space,  stability,  convenience,  and  illuminated  area, 
most  laboratory  wave-tank  scattering  measurements  are  made  at  ranges  R0  within  the 
Fresnel  region.  This  region  extends  from  several  wavelengths  "rom  the  antenna  out  to 
**=>2 D2/kp,  where  D is  the  antenna  diameter.  The  Fresnel  region  of  a focussed  parabolic 
antenna  is  particularly  suited  to  laboratory  scattering  measurements.  The  microwave 
antenna  pattern  at  Rq  will  be  derived  not  for  exact  computation  but  to  show  its  similar- 
ity to  any  other  focussing  system  and  to  obtain  important  scaling  parameters.  Refer  to 
Fig.  1 for  the  coordinate  geometry. 

Following  Silver  [5] , the  microwave  field  at  a point  (,r,  y)  on  the  water  surface 
can  be  written  (suppressing  an  exp  (/cut)  time  dependence)  as 
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£\x,  7)  = ~ j fXu,  v)  (r.  + r2rl  erp  [-Oto^i  + r2)J  dudv.  (9) 

Here  F{u,  v)  is  the  antenna-feed  field  distribution  over  the  aperture  A.  Expanding  the 
phase  term  to  second  order  in  X.R51  sin  0,  yRo1,  uRp~l,  and  vRp*  gives 

ri  + r2  * Ro  * f + e + x cos  6 + ~ (u2  + v2  )L(e) 

mt 

- uxJio1  sin  0 vyR  $ 

+ x2(2Rp)1  sin2  6 + y2(2Rp)^ , (10) 

whe»3  higher  order  terms  have  been  neglected.  The  coefficient  of  the  u2  + v2  term  is 


+ 


1 

f+e 


(ID 


where  c is  the  spacing  between  the  antenrs  feed  and  the  “infinity”  focal  point  of  the 
parabola,  as  shown  in  Fig.  1.  A focussed  condition  is  defined  by  an  e satisfying  the  lens 
equation 


Ue)  = 0.  (12) 

To  this  order  of  approximation,  phase  variations  at  (x,  y)  = (0,  0)  due  to  a finite  aper- 
ture diameter  can  be  “nulled”  by  a proper  choice  of  e.  The  parabola  can  be  effectively 
focussed  on  the  water  surface  and  can  produce  approximately  plane  waves  over  the  illumi- 
nated area,  even  at  short  range. 


Using  Eq.  (10)  in  the  phase  term  of  Eq.  (9)  and  approximating  the  denominator  by 
(l*o  + f)  gives 


£(x,  y)  = 


ik0 

2*IR0  * f ) 


e*P  jn*o  («o  + f + 


+ x cosfl  + 


c2  sin20 
2R0 


* f '*%*)]* 

x J_  „p  [-»o(f;  * ^)]  ■“> 


The  feed  gain  function  F(u,  v)  has  been  written  as  F\(u) Fplv)  and  defined  as  F\  - 
Fz  - 0 outside  of  the  antenna  aperture,  which  is  adequate  for  this  discussion.  Both 
integrals  in  Eq.  (13)  are  a Fourier  transform  of  a product  of  two  functions,  that  is,  a 
convolution  of  two  Fourier  transforms.  The  coordinate  transform  pairs  are  u,  ax  sin  6 
and  v,  cry,  where  a = kp/Rp.  Denoting  a Fourier  transform  by  . and  a convolution  by  *, 
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Eq.  (13)  becomes 


£(*.>  > = Co  sxp{-ikox  co»B)fl(x)f2(y), 


(14) 


where 


Co  » 


ik  o 


exp  [i*o(^o  + f + f)l. 


2x(i?o  + f) 

f\(x)  = exp  ^-i^x2  sin2  J^1(ttr  an  8)  * exp 

*(- 


(15) 


rtcr2x2  sin2  d 


and 


fiiy)  = exp 


koL(e) 


(•‘I ■ yt)  [ 


)]■ 


( itt2y2  \ 
F*{ay)  * ”P * 


(16) 


(17) 


Exact  calculation  of  the  fields  * tedioic.  and  requires  the  feed  gain  function  to  be 
specified.  The  direct-measurement  techniques  discussed  in  section  4 are  more  practical  and 
more  accurate.  Th?  important  qualitative  feature  is  that  for  L(e)  ~ 0,  the  antenna  pattern 
E(x,  y)  becomes  approximately  the  Fourier  transform  of  the  feed  gain  function,  scaled  by 
a sin  6 or  a.  Estimates  of  die  minimum  real-space  pattern  widths,  in  both  direction,  are 


Ax  *= 


1 ^ E0Xq 

a sin  0 Au  2rD  sin  0 


(18) 


and 


Ay  ' 


aAv 


Ro^o 

2 wD 


(19) 


Here  Au  and  Au  are  the  widths  of  Ft  (u)  and  F 2(0)  respectively  and  are  approximated  by 
D In  practice,  most  feeds  have  a gain  function  tapered  from  the  center  to  edge  of  the 
parabola.  This  reduces  An  ar.d  Ac  but  improves  the  sideiobes.  Other  choices  of  e further 
broaden  the  real-space  pattern.  The  pattern  width  in  each  direction  in  this  case  can  be 
estimated  using 


and 


(20) 


(21) 
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From  another  viewpoint,  the  radius  of  the  first  Fresnel  zone  centered  on  the  "spot” 
illuminated  at  the  water  surface  is  given  by 


rF 


(22) 


Taking  the  ratio 

* V*Tq 

rF  Au  D 

reveals  that  for  sufficiently  short  ranges  this  first  Fresnel  zone  is  under -illuminated.  In 
this  sense  one  has  an  incident  plane  wave,  and  the  illuminating  patterns  are  transforms,  as 
discussed  previously. 
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3.  PREDICTED  PROBE  RESPONSE 


The  successful  use  of  this  wave  probe  depends  on  understanding  the  effects  of  the 
illumination  pattern  on  the  beckscattered  poorer.  Besides  illuminated  area,  important 
factors  are  resolution  in  water  wavenumber  space,  and  absolute  calibration  in  terms  of 
wave  height.  For  cases  when  the  an  term:  i pattern  is  important,  the  backscattered  power 
is  best  calculated  using  an  iterative  integral  method  [6]  which  is  particularly  simple  in  first 
order.  This  is  done  in  detail  in  the  appendix.  The  fnst-o»der  retiimed  power  is  proportional, 
as  described  in  the  appendix,  to  a dimensionless  cross  rection,  or  cross  section  per  unit 
illuminated  area: 


o0  - fc(0)l2 ♦<**,*,, 


" = <*B,  0) 


The  transfer  coefficient  lg(0)!2  depends  on  the  polarization  and  on  the  incidence  angle. 

It  is  given  in  the  appendix  and  is  plotted  by  Wright  [7].  The  water  wave  spectrum  to  be 
measured,  ’J f(feg,  0,  w),  is  convolved  with  an  antenna  response  function  V(feg,  0)  which  is 
normalized: 


V(kx,ky)dkxdky  = 1.  (25) 

This  function  determines  the  resolution,  Akx  and  Aky,  possible  about  k = (kg,  0).  An 
interesting  feature  u>  that  V(kx , ky ) is  related  to  the  spatial  Fourier  transform  of  the 
antenna  pattern  at  the  water  surface  and  thus  to  the  feed  gain  function.  More  explicitly, 

v{kx,ky)  = ’}(feJC,ky)*T7(fcx,ky),  26) 


where 


a sin  & 


(27) 
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Tots  inverse-transform  relationship  between  fetd  gain  function,  antenna  pattern,  and 
fc-space  response  is  typical  of  focussim  systems.  The  scale  factors  in  Eq,  (27)  allow 
one  to  predict  resolutions  for  other  sets  of  parameters  D , 0 k$,  and  Rq,  given  one 
result  initially. 

The  choice  of  antenna  and  range  usually  involves  a compromise.  If  a large  tntenna 
diameter  and  short  range  are  chosen  to  give  severe  underillumination  and  a good  plane 
wave,  the  focussed  spot  may  be  only  a Bragg  wavelength  or  less  in  diameter  and  yield 
pocr  resolution  in  wavenumber.  If  overillumination  is  chosen,  the  Fresnel  terms  will 
produce  sidelobes  in  V(kx.  ky ) and  broaden  the  main  lobe.  For  the  antennas  we  have 
investigated,  maximum  resolution  in  wavenumber  occurs  near  the  L(e)  = 0 condition. 

It  is  practical,  however,  to  enlarge  the  illuminated  area  by  limited  defocussing. 

In  any  case,  exact  calculation  of  V(fex,  fey)  is  complex  and  depends  critically  on  the 
phase  as  well  as  on  the  amplitude  of  the  fields  produced  by  the  antenna  feed.  It  is  much 
simpler  and  more  accurate  to  measure  V(fex,  ky)  directly  using  monochromatic,  mechan- 
ically generated  water  waves.  In  this  caoe  the  wave  spectrum  is  narrow  compared  to 
V(kx,  ky).  If  the  backscatte;ed  power  is  measured  while  the  wave  frequency  (or  kx>  is 
slowly  scanned,  F[fe(cn),  0}  is  traced  out.  Although  the  cross  section  is  not  independent 
of  range  and  antenna  properties,  it  is  given  by 

a(u>)  - &irk04h2  i*2(0\ lv[fe(W),OMo,  (28) 

where  h h the  mechanical-wove  rms  amplitude  and  A a is  tire  illuminated  area,  as  denned 
in  the  appendix.  To  make  this  measurement,  it  is  only  necessary  to  keep  h constant;  the 
magnitude  of  V is  obtained  from  the  normalization  condition,  Eq.  (25). 


4.  PROBE  RESPONSE  MEASUREMENTS 

The  fe-space  response  of  this  probe,  or  V[«(w),  0] , is  measured  by  a technique 
which  we  believe  has  not  been  previously  reported.  With  an  antenna  mounted  above 
the  water  as  in  Fig.  1,  monochromatic,  mechanically  generated  water  waves  are  propa- 
gated in  the  x direction  through  the  illuminated  area.  The  baokscattered  power,  doppler 
shifted  at  the  water  wave  frequency  , is  plotted  as  a function  of  cjw,  which  is  linearly 
swept.  A sketch  of  the  loudspeaker-driven  mechanical  wavemaker  is  shown  in  Fig.  5.  A 
linear  voltage  ramp  drives  a voltage-controlled  audio  oscillator  whose  output  is  fed  to  the 
speaker  after  amplification.  A tew  experimental  pitfalls  should  be  mentioned.  The  dimen- 
sion b in  Fig.  5 must  be  chosen  to  avoid  mechanical  resonances  in  the  driven  water  col- 
umn, especially  near  a Bragg  frequency,  An  excessive  insertion  depth  a severely  reduces 
amplitude  output.  At  frequencies  nea”  8 Hz,  amplitude  fluctuations  are  produced  which 
vary  slowly  over  several  minutes.  These  are  perhaps  due  to  a resons'  wave-wave  inter- 
action mechanism  (S] . 

A representative  trace  of  backscattered  power  vs  we. tor  wave  frequency  fy  is  shown 
in  Fig.  6.  For  this  example,  these  parameters  apply:  VV  polarization,  D = 30.5  cm, 

R0  ~ 1.22  m,  = 9.375  GHz,  6 - 35°,  and  \g  ~ 1.95  cm.  The  pattern  shape  is  near 
gaussian,  or  parabolic  when  plotted  logarithmically  (Fig.  6).  This  gives  very  steep  filter 
skirts  indeed  and  makes  the  teehnicue  useful  as  a wave  probe.  Frequency  resolution 
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can  be  read  directly  (for  example,  using  a full  width  at  6 dB  below  maximum 
criterion).  Wavenumber  resolution  is  obtained  from 


Akx  = AcJjy 


dco(fefl)  -1 

3fe  J ’ 


(29) 


where  d<o(fes)/dfe  is  evaluated  using  Eq.  (8).  Structure  appearing  in  the  trace  at  f < fg 
is  probably  caused  by  nonsinusoidal  waveforms.  A low-frequency  wave  could  have  a 
spatial  harmonic  at  k&  ■ This  trace  is  one  of  a series  taken  under  identical  conditions, 
except  that  B was  varied  in  5°  increments  for  each  trace.  As  discussed  in  section  1,  the 
measured  peak  values  of  are  plotted  in  Fig.  3 vs  the  Bragg  resonant  frequencies 

predicted  using  Eqs.  (1)  and  (8).  The  peak  value  of  V(kx  ) should  be  inversely  propor- 
tional to  the  illuminated  area,  which  in  turn  varies  as  sin  9.  This  is  verifed  in  Fig.  7 where 
the  integrated  area  under  the  main  peak  of  each  trace  is  plotted  vs  sin  6.  A straight-line 
dependence  is  verified.  As  6 approaches  90°,  interference  from  the  large  specularly  re- 
flected beam  becomes  severe.  Thus  the  data  have  more  scatter  for  these  angles. 


Fig.  7— Aren  under  a aeries  of  antenna  response  curve*  like  the 
one  in  Fig.  6,  vs  ain  9.  Resolution  scales  inversely  with  the 
length  of  the  illuminated  area  and  directly  with  sin  6. 


By  taking  traces  at  various  antenna-feed  positions,  the  dependence  of  fe -space  resolu- 
tion on  e can  be  investigated.  Results  of  two  series  of  traces  are  given  in  Fig.  3 and  9. 
Parameters  for  Fig.  8 are  fo  = 4.30  GHz,  D = 61  cm,  6 = 60°,  = 7.00  cm,  and  a focal 

length  of  20.3  cm;  for  Fig.  9,  f0  = 9.375  GHz,  D = 30.5  cm,  9 = 45°,  XB  = 2.26  cm,  and 
a focal  length  of  10.2  cm.  Both  cases  had  the  same  range  Rq  of  1.22  m.  For  both 
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FEEO  POSITION  (cm) 

Fig.  8— Resolution  Afcx  vs  feed  position  (arbitrary  origin)  for  0 
= 60°,  D = 61  cm,  = 7.00  cm,  R<>  * 1.2  m,  and  a focal  length 
of  20.3  cm. 


FEED  POSITiONIcm' 

Fig.  9— Resolution  vs  feed  position  (arbitrary  origin)  for 
6 * 45°,  D - 30.5  cm,  \g  * 2.26  cm,  Rq  ~ 1.2  m,  and  a focal 
length  of  10.2  cm. 
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figures  the  wavenumber  resolution  obtained  from  Eq.  (29)  is  plotted  vs  feed  position  in 
cm.  The  criterion  for  Ate#,  is  full  width  at  6 dB  below  maximum.  The  abscissa  in  each 
case  has  an  arbitrary  origin  because  for  each  feed  the  effective  ray  origin  position  is  not 
known.  The  polarizations  used  arc  indicated.  In  both  cases  a range  of  ±15%  in  e about 
the  focussed  position  is  available  for  varying  the  illuminated  area,  without  severely  degrad- 
ing the  resolution.  The  best  wavenumber  resolution  from  the  data  given  here  is  Akx/k^ 

= 10%  for  Fig.  9 and  Akx/k g = 20%  for  Fig.  8. 

Several  real-space  antenna  patterns  were  measured  to  check  their  dependence  on  e 
and  Rq  and  to  determine  the  illuminated  area  /to.  defined  by  Eq.  (46).  These  two-way 
patterns  are  measured  by  vertically  lowering  a metal  sphere  of  known  cross  section  [9a] 
through  the  radiated  antenna  field,  with  the  axis  of  symmetry  of  the  antenna  horizontal. 
The  range  Ro  was  1.22  m.  The  diameter  of  the  metal  sphere  was  typically  ~ Xo-  Care 
in  boresighting,  or  aligning  the  sphere  to  pass  through  the  center  of  the  pattern,  is  neces- 
sary since  the  main  beam  direction  can  deviate  from  the  apparent  antenna  axis  of  sym 
metry  by  a few  degrees.  Pattern  widths  in  the  E and  H planes  are  plotted  vs  feed  position 
in  Fig.  10  for  the  30.5-cm  parabola  with  focal  length  10.2  cm  at  9.375  GHz.  Thus  each 
plotted  point  is  obtained  from  a measured  pattern  whose  width  is  taken  as  the  full  pattern 
width  at  6 dB  below  maximum.  Figure  11  gives  the  same  measurements  for  a 61 -cm 
parabola  with  focal  length  20.3  cm  at  4.30  GHz.  In  both  cases  the  shape  of  the 
plotted  data  has  a broad  minimum  defining  the  focussed  position.  Because  of  the  near- 
gaussian  patterns,  the  illuminated  area  is  proportioned  to  the  product  of  these  two  widths. 
As  a check  on  depth  of  field.  Fig.  12  shows  the  E-plane  main-lobe  width  vs  Rq  for  the 
4.30  GHz  case.  The  feed  position  e was  chosen  to  give  a focus  at  a range  of  1.5  m in 
this  case.  The  pattern  width  scales  linearly  with  Rq,  as  predicted  by  Eq.  (18).  The 
depth  of  fie  ri  is  thus  adequate  for  most  water  wave  measurements  made  in  a laboratory 
tank.  Illuminated  area  at  any  0 can  be  obtained  from  the  measured  90°  values  by  divid- 
ing by  sin  8.  The  pattern  did  not  deteriorate  drastically  either  in  main-lobe  shape  or  in 
sidelobe  level  for  at  least  a ±15%  range  in  e//  about  the  focus  for  all  frequencies  and 
parabolas  used.  Thus  changing  illuminated  area  by  focussing  a parabola  in  the  near  field 
is  a viable  procedure.  By  comparing  real-shape  pattern  widths  for  the  E and  H planes, 
the  resolution  A ky  can  be  obtained  from  Akx. 

Finally,  as  a consistency  check  on  all  of  the  calibration  measurements  made  using 
this  probe,  returned  power  scattered  from  a monochromatic  mechanically  generated  wave 
at  k = kjj  was  measured  absolutely  at  several  angles;  these  measurements  are  compared  in 
the  appendix  with  the  results  of  first-order  scattering  theory.  The  data  are  taken 
from  the  same  series  of  9.375-GHz  traces  taken  at  5°  increments  in  9 previously  discussed. 
The  quantities  Aq,  V(ks,  0),  and  h were  measured  independently  for  each  angle  and  used 
in  Eqs.  (28),  (43),  (44),  and  (45)  to  predict  PT/Pt  vs  9.  These  predicted  values  are  com- 
pared to  actual  measured  values  of  Pr/Pt  in  Fig.  13,  where  both  are  plotted  vs  6.  The 
measured  returned  power  is  compared  in  each  case  to  a standard  reflection.  This  in  turn 
is  calibrated  with  the  return  from  a metal  sphere  at  a known  range  and  focus  setting.  The 
measured  fc-space  patterns,  together  with  the  normalization  condition,  Eq.  (25),  give  the 
V(ku,  0)  term.  The  water  wave  amplitude  h was  measured  in  each  case  by  photographing 
the  image  of  an  illuminated  slit  reflected  in  the  water  surface.  Both  measured  and  calcu- 
lated values  of  Pr!Pt  «*e  normalized  to  a water  wave  amplitude  of  0.27  mm,  peak  to  peak. 
Actual  amplitude  was  kept  within  a factor  of  two  of  this  value  for  each  case.  The  value  of 

6 = ksh  sin  9 = 2kGh  sin  6 (30) 
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Pig.  11 — Rfil-tptH  antenna  pattern  width  vs.  feed  position 
(same  origin  as  Fig.  8).  D * 61  cm,  foca)  length  - 20.3  cm,  f0 
* 4.30  GHz,  and  Rq  ' 1.2  m.  Both  polarization  directions  are 
given,  as  indicated. 
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Fig.  1 2 — Pattern  width  vs  range  R0  for  a fixed  feed  position. 
D - 61  cm.  focal  length  = 20.3  cm,  and  /0  ■ 4.30  GHz. 


Fig.  13— Measured  value*  of  backacattered  power  vs  0.  The 
■olid  line  give*  the  values  predicted  by  the  first-order  scattering 
theory  in  the  appendix. 
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has  a maximum  of  Q.G6.  Thus  the  Rayleigh  criterion,  6 < 1,  was  certrinly  satisfied,  and 
first-order  scattering  theory  should  be  sufficiently  accurate.  Water  wave  attenuation  is 
measured  directly  by  using  two  illuminated  slits.  If  care  is  taken  to  blow  off  contaminant 
surface  films,  attenuation  is  small  enough  to  be  neglected  over  the  length  of  the  illuminated 
area.  This  area  /to  is  obtained  from  real-space  patterns,  as  described  previously.  The 
transfer  coefficient  lg(0)l2  is  calculated  using  55(1  - i0.55)  for  the  relative  dielectric 
constant  of  water.  Agreement  is  satisfactory.  At  the  angles  nearer  normal  incidence, 
interference  from  the  specularly  scattered  beam  probably  caused  the  measured  values  to 
be  slightly  high.  However,  the  results  of  first-order  scattering  theory  are  certainly  bom 
out,  and  the  calibration  procedures  described  in  this  report  are  consistent. 
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APPENDIX 


CALCULATION  OF  BACKSCATTERED  POWER 

In  this  case  where  the  antenna  pattern  is  important,  the  backscattered  power  is  best 
calculated  using  an  iterative  integral  method  which  is  particularly  simple  in  first  order  [7] 
Referring  to  Fig.  1,  the  resultant  total  field  in  all  space  is  denoted  by  Sq.  Hq  an  un- 
perturbed water  surface.  When  the  surface  is  now  perturbed  by  a displacement  'fix,  y.  f), 

ossume  that  the  resultant  total  field  in  all  space  can  be  written: 

E = E0  + «i  + *2  * (31a) 

H = Hq  + hj  + hj  + ...  . (31b) 


Successive  terms  are  proportional  to  ascending  powers  of  a dimensionless  expansion  pa- 
rameter proportional  to  surface  displacement,  for  example,  the  5 defined  in  Eq.  (30). 
Using  a genera,'  scattering  theorem  (9b) , the  average  power  backscattered  by  the  rough 
surface  y(x,  y,  t)  is 


Pr 


1 

16P,2 


r 


E X Hq  - Eq  X H) 


n dc'j2. 


(32) 


where  P,  is  the  transmitted  power.  Here  n is  ho  local  surface  normal,  and  s'  is  a surface 
of  integration  enclosing  y(x,  y,  f)  above  and  below,  as  shown  in  Fig.  1.  Kerr’s  method 
is  readily  adapted  to  calculation  of  the  received  power  spzctral  density  or  doppler  spec- 
trum mu j) 
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jinx  E)  * AH0  + (nX  H)  • AEq]  ds  , 


(33) 


where  AE0,  A‘ln  are  the  discontinuities  in  the  unperturbed  field  across  the  perturbed 
surface  y(x,  y,  t)  and  n is  the  normal  to  the  perturbed  surface. 


AE0  = (E0**r  - E0w*ler) 

* = o 
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and  timiiariy  for  AHo  . Neglected  in  Eq.  (33)  were  terms  of  order  (Au/uii'  where 
Acj  is  the  width  of  the  doppler  spectrum  and  wo  ® the  microwave  frequency.  The 
finite  antenna  pattern  assures  negligible  contribution  to  the  spatial  integral  at  large 
distances.  The  total  backscattetcri  power  is 

Pr  ~ f D(cj)duj.  (35) 


The  zero-order  fields  at  the  surface  are  obtained  in  principle  from  Eq.  (14).  However, 
these  fields  can  be  assumed  locally  plane  for  the  purpose  of  computing  AEq  md  AH0. 
That  is,  the  relationships  between  electric  and  magnetic  fields  and  their  derivatives  are 
assumed  to  be  those  for  plane  waves.  This  may  introduce  errors  in  the  received  power 
of  magnitude  (D/R^  )2. 

Substituting  Eq.  (34)  in  Eq.  (33)  reveals  the  advantage  of  this  calculation  technique: 
•i , h j , and  hitter  order  terms  are  not  needed  to  evaluate  the  integral  to  first  order. 
Equation  (33)  becomes 


Z?fa>) 

P, 


X 


dxdy 


(36) 


where  kg  is  2*o  cos  6.  Tht  product  £,( x,  y)W,(x,  y)/2  is  the  complex  illumination  pat- 
tern, where  Hi  = Ej/Zo,  with  Z0  the  characteristic  impedance  of  free  space.  Thus 
i£,(0,  0)  H,i 0,  0)1/2  is  the  mean  incident  intensity  at  the  center  of  the  illuminated  area. 
The  angular  dependence  is  giver,  bv 

2 

vv 

and 

fffj  = 

where  f ’ is  the  dielectric  constant  of  water.  The  subscripts  define  the  E-field  polarization; 
for  example,  VV  impLcs  polarization  which  is  vertical,  or  in  the  plane  of  incidence,  for 
both  transmit  and  receive.  Note  that  in  actual  CW  doppler  measurements,  an  audio- 
frequency offset  is  introduced  to  prevent  folding  positive  md  negative  velocities  onto 
one  another  in  the  detection  process.  This  offset  has  been  omitted  here. 

As  in  mo6t  scattering  problems,  it  is  useful  to  introduce  the  scatteriiig  cross  section 
o,  which  is  defined  by  the  radar  range  equation 
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Again,  by  definition  the  antenna  gain  G is 


G - 4**02  »» 


2P, 


(40) 


and  Eq.  (40)  can  be  rewritten 


T,  ’ V ^ 


(41) 


The  definition  of  o based  on  Eq.  (39)  is  valid  in  the  far-field  limit  Rq  > D2!\q . It  i£ 
only  then  (and  in  the  geo  metrical -optics  limit)  that  a is  truly  independent  of  range  and 
antenna  properties.  For  our  case,  however,  we  substitute  Eq.  (4i)  in  Eq.  (36)  and  replace 
E(R0)H(R0)  by  E,(0,  0)/f(0,  0).  This  gives  a definition  of  a applicable  for  this  near-held, 
doppler-radar  case: 


o(cj)  ~ lim 
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(42) 


where 


a 


ZO(h>)  dui. 
» 


(43) 


Tb  definition  of  a is  extremely  useful  because  the  simplest  way  to  calibrate  received 
power  in  these  measurements  is  to  compose  it  with  the  btctocstoed  power  from  a target 
of  theoretically  known  cross  section,  usually  a small  metal  sphere.  The  theoretical  cross 
section  is  computed  bom  a boundary-value  problem,  assuming  incident  plane  waves.  As 
this  is  consistent  with  our  assumptions,  no  further  approximation  is  made  when  using  this 
calibration  procedure.  To  summarize,  by  comparing  the  backseat  tried  power  from  the 
rough  surface  with  that  from  a highly  conducting  sphere  placed  at  the  center  of  illumine 
tion,  we  can  obtain  a crow  section  (in  cm2,  say)  for  the  rough  surface. 

If  y(x,  y.  0 is  * random  surface  for  which  a suitable  ensemble  can  be  constructed, 
a dimensionless  cross  section  o0  can  be  defined: 


®0 


(44) 


where 
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A0  = 


EjfayWite  y) 
E,i 0,  0)Hi(Q,  0) 


2 

dxdy 


(45) 


Thus  oo  “ * dimensionless  cross  section  normalized  by  the  illuminated  area  defined  by 
Eq.  (45)  and  is  independent  of  range  and  antenna  properties  in  many  practical  cases.  To 
compute  Oq,  we  use  the  identity 

I///— H’ s fJJfJfo*- v- r' 

y f{x  + Ax,y  + Ay,  t * t) 

X dxdy<Ud\xd\ydT  (16) 

and  take  an  ensemble  average  of  Eq.  (42).  The  result  is 

a0(co)  - W(«),2Jj(jf^Ax,A,,r)ir(Ax.  Ay) 


X e**aA,  e*<*x  d\x  d\y  dt , 

(47) 

where  Z( Ax, 

Ay,  r)  is  the  surface  covariance 

Z{ Ax,  Ay.r)  = (y(x,  y . f)7(x  + Ax,y  + A,,  1 ♦ r>). 

(48) 

*(AX,  A,)  = A jfjfW  y ) V*{x  * Ax,  y + Ay)  dxdy. 

09) 

and 

mi  - **<*•  y> 

' ’ E.to,  0)«,<0,  0) 

(50) 

Now  with 

+(*x,  *„  w)  = £ 

(51) 

and 

V(*a,  *,)  = £ 

(52) 

we  can  rewrite  Eq.  (47)  as 

O0<«)  = 16»kf)4  lg(0)l*  *(**,*,,  <•>)*  V(kx,  ky)j 

1*  m (*a.  0) 

(53) 

which  is  diacumed  in  section  3. 
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